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ABSTRACT 



Herschel-KWl observations of water in the low-mass star-forming object L1448-MM, known for its prominent outflow, are presented, 
as obtained within the 'Water in star-forming regions with HerscheV (WISH) key programme. Six H'*0 lines are targeted and detected 
(£up/ytB -50-250 K), as is CO J= 10-9 {E^^lk^~2Q5 K), and tentatively O lio-loi at 548 GHz. All hnes show strong emission 
in the "bullets" at |[)|>50 km s"' from the source velocity, in addition to a broad, central component and narrow absorption. The 
bullets are seen much more prominently in HtO than in CO with respect to the central component, and show little variation with 
excitation in H2O profile shape. Excitation conditions in the bullets derived from CO lines imply a temperature >150 K and density 
> 10' cm"\ similar to that of the broad component. The H2O/CO abundance ratio is similar in the "bullets" and the broad component, 
-0.05-1.0, in spite of their different origins in the molecular jet and the interaction between the outflow and the envelope. The high 
H2O abundance indicates that the bullets are H2 rich. The H2O cooling in the "bullets" and the broad component is similar and higher 
than the CO cooling in the same components. These data illustrate the power of Herschel-WFl to disentangle different dynamical 
components in low-mass star-forming objects and determine their excitation and chemical conditions. 

Key words. Astrochemistry — Stars: formation — ISM: molecules — ISM: jets and outflows — ISM: individual objects: L1448 



1. Introduction 

Low-mass star formation is accompanied by the launch of pow- 
erful jets that dri ve strong shocks into the parental material (e.g., 
lArce et aT ] l2007h . both in the form of shell shocks that interact 
with the molecular envelope, and in the form of the jet itself 
(e.g., Inirano et al. 2010). The jets are generally not uniform, 
but consist of small condensations ("bullets'Q), that can be at- 
tributed to internal working surfa ces caused by episodic ejec - 
tion of material from the protostar (ISantiago-Garcfa et al .1120091) . 
The bullets are chemically rich, in particular in oxygen-bearing 
species (Tafalla et al. 2010), and should therefore display high 
abundances of water. H2O is one of the best shock tracers be- 
cause the abundance can increase by orders of magnitude up 
to <10 through sputtering of gra in mantles and formation in 
the g as phase at high temperatures ("Flower & Pineau des ForetsI 
I2OIQ) . Also, due to the large velocity gradients in shocks, emis- 
sion is effectively o ptically thinner than in a quiescent object 
dFranklin et al.ll200 8). 

H2O has previously been observed in outflowing gas from 
young stellar objects (YSOs) using, e.g., Odin, SWAS, and 
ISO. The Heterodyne Instrument for the Far-Infrared (HIFI) on 
Herschel opens up spectrally resolved, high angular resolution 



obser vations of both H7O a nd high-/ CO in low-mass protostars 
(e.g. lLeflochetal.1 1201 Ol lYildiz et all [20 Id iKristensen et"al] 
l2010h . The latter paper presents H2O observations of the low- 
mass star-forming region NGC1333. The H2O line profiles are 
complex, consisting of three velocity components characterised 
by their width. The broadest component (At;>20 km s"') orig- 
inates in the interaction between the outflow and the envelope; 
the medium-broad component (5<Ay<20 km s"') arises where 
the jet impacts on the inner, dense envelope; a narrow absorp- 
tion component (At;<5 km s"') is attributed to the cold outer 
envelope. These initial HIFI results show that H2O is one of the 
best dynamics tracers in YSOs, and confirm that the H2O/CO 
abundance increases with velocity to ~ 1 . To date, H2O has not 
been uniquely identified in the bullets, and its abundance with 
respect to other jet species is unknown. Moreover, the different 
nature of the shocks could lead to different H2O abundances. 



L1448-MM is a low-mass 
in Perseus (L=11.6 L q, Menv 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 

' Bullets and extremely high-velocity (EHV) gas are used inter- 
changeably throughout. 



Class protostar located 
=1.9 Mo, d = 235 pc; 
Ivan Dishoeck et all |201 1 *). This object is the prototype of the 
class of sources showing bullet emission in CO (Bachille r et al.l 
Il990h . with red- and blue-shifted emission peaks at ~ ±50 km 
s'^ with respect to dlsr = +5.2 km s"'. The bullet s are clearly 
visibl e in other shock tracers such as SiO (e.g. iNisini et aLl 
2007), and have been imaged at high spatial resolution using 
interferometers (e.g., Guilloteau et al. 1992, Hirano et al. 20 1_^. 
Our observations of H2O within the framew ork of the 'Water in 
star-forming regions with HerscheV (WISH: Ivan Dishoeck et alJ 
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1201 Ih key programme further explore the chemical and excita- 
tion conditions of shocked gas in low-mass star-forming regions. 

2. Observations and results 

The central position of L1448-MM (03''25"'38 ^9; -h30°44'05'.'4; 
J2000) was observed with HIFI on Herschel dde Graauw et al.l 
l2010llPilbrattetal.ll2010h in seven different settings covering 
six H^^O, two H^^O and one CO transition {E^/kB - 50- 300 
K; Table [3] available online). Data were obtained using the dual 
beam-switch mode with a nod of 3', except for the ground-state 
ortho-H20 line at 557 GHz, where a position switch to (+10'; 
+5') was used. The diffraction-limited beam size ranges from 
19" to 39" (4500-9500 AU). Data were reduced using HIPE ver 
4.0. The calibration uncertainty is taken to be 10% for lines ob- 
served in Bands 1, 2, and 5 while it is 30% in Band 4. The point- 
ing accuracy is ~2". A main-beam efficiency of 0.65-0.75 is 
adopted (Table |3]l. Subsequent analysis of the data is performed 
in CLASS including subtraction of linear baselines. H- and V- 
polarizations are co-added after inspection; no significant differ- 
ences are found between the two data sets. '^CO 3-2 and '"^CO 
3-2 were observed with the James Clerk Maxw ell Telescope 
(JCMT). CO 2-1 was observed a t the IRAM 30m (iTafalla et alJ 
|20T0 ). and CO 4-3 at the JCMT (iNisini et alJIiOOOl) . 

All targeted Unes are detected, except Hj^O Ih-Ooo, ob- 
served with the main isotopologue (Fig. [TJ. This figure strik- 
ingly demonstrates that the bullets are much more prominent 
in H2O than in CO relative to the central component. Also, the 
H2O profiles are very different from those seen in NGC 1333 
dKristensen et al.i r2010). Each line consists of multiple compo- 
nents: two Gaussian peaks at velocities of ~ ±50 km s"' with re- 
spect to Usource; a broad (~50 km s"') component centered close 
to Usource; and a narrow (~5 km s"') component seen in absorp- 
tion at ^source- The extremely high-y emission (henceforth EHV-B 
and EHV-R) is detected in all lines, although the detection in the 
li\^0 lio-loi line is tentative (2cr; 0.13 K km s"'). The emis- 
sion components are well reproduced by Gaussians and fitted 
to obtain the integrated intensity (Table [TJ. The H2O line pro- 
files vary little with excitation. The shape of the EHV component 
in the CO profiles resembles that of H2O (Fig. lA.ll appendix), 
whereas the the broad component is much weaker in CO. The 
CO lines also contain a narrower component (~10 km s ') not 
seen in H2O. 

To compare observations done with different beam sizes, all 
observations are scaled to a common beam of 22", the beam at 
98 8 GHz. The s caling follows the recipe outlined in Appendix B 
of iTafalla et al ](|2010), where a linear scaling with beam size is 
found to be appropriate for the EHV components and a power- 
law scaling with index 0.75 for the broad component. It is only 
significant for the 557 GHz lines due to the larger beam (39"). 

The tentative detection of the VL\^0 lio-loi (548 GHz) fine 
is used to constrain the optical depth, t, of the main isotopo- 
logue, assuming that the line itself is optically thin and that 
['*'O]/[''^O]=550. The optical depth, t, is -25 in each EHV com- 
ponent and 9 in the broad component implying that emission 
from Hj^O 1 lo-loi is optically thick. The 3cr upper limit on the 
H^'^O lii-Ooo fine impHes t<10 and <30 in H^''0 for the broad 
and EHV components, respectively. '^CO 3-2 emission is not 
detected in any of the components. The '^CO 3-2 optical depth 
is <2 and <0.3 in the EHV and broad components, respectively, 
implying that the CO emission is optically thin. 

The integrated intensities of each of the EHV and broad com- 
ponents are used to construct Boltzmann diagrams, where the 
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Fig. 1. Continuum-subtracted HIFI H2O, R^^O and CO spectra 
obtained at the central position of L1448-MM (usomce - 5.2 km 

Table 1. H2O, ii\^0 and CO emission in each velocity 
component" . 



EHV-B EHV-R Broad 
Transition rms'' j Tmb Tmb j Tub dti 
(mK) (Kkms-') (Kkms'') (Kltms'') 



H2O lio-loi 


9 


2.54 


3.00 


14.4 


lii-Ooo 


11 


3.98 


4.00 


14.8 


2o2-lii 


14 


4.01 


3.86 


17.6 


2ii-2o2 


14 


2.47 


2.53 


9.9 


3l2-3o3 


13 


2.34 


1.14 


7.8 


3 12-221 


66 


6.22 


3.32 


7.9 


Hj*0 lio-loi' 


3 


0.13 


0.11 


0.23 


lii-Ooo'' 


11 


<0.25 


<0.I9 


<0.28 


CO 2-1 


51 


2.79 


4.12 


38.5 


3-2 


39 


4.82 


6.58 


44.9 


4-3 


160 


11.9 


15.1 


74.6 


10-9 


66 


7.30 


8.42 


17.8 


EWKM' (kms-') 




38.4±3.2 


22.4±3.0 


47.6±4.5 


ulsr' (kms"') 




-43.4±1.3 


59.3±1.8 


16.8±2.2 



Notes. Obtained from Gaussian fits to each component. All Gauss- 
fit parameters are provided in Table [2] available online. Measured 
in 1 kms"' bins. '•'^^ Obtained by fixing the position and width of 
Gaussian functions and only fitting the intensity. Upper limits are 
3cr. Intensity- weighted average of values determined from Gaussian 
fits to the H2O lines. Uncertainties include statistical errors only. 



column density of each upper level in the optically thin limit 
(A^*p) per sub-level is plotted versus E^p (Fig.|2]i. In the diagrams, 
A^*p for the H2O ground-state transitions is lower than expected, 
confirming that the true A^up is higher and that the lines are op- 
tically thick. Excluding these two lines from fuither analysis, 
the resulting rotational temperature is remarkably similar for all 
three components (~45-60 K), while the beam-averaged column 
density of the broad component is ~5 times higher than that of 
the EHV components (~ 10 and 2x10'^ cm"^, respectively; Table 
|2]l. Rotational diagrams are also made for CO and results are 
given in Tabled T^ot is significantly higher, 150 K in the EHV 
components and 80 K in the broad. The column density is also 
higher, -10"' cm -, and the inferred H2O/CO ratio is -0.002. 
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Fig. 2. Rotational diagrams of the EHV and broad components 
for H2O (red circles) and CO (blue triangles). CO has been 
shifted downwards by a factor of 100. The ground-state tran- 
sitions (grey) are excluded from the fit (dashed lines). Hl^^O in- 
tensities x550 are shown as open circles. 

Table 2. Rotational-diagram results for each velocity compo- 
nent. 







H2O 




CO 


Component 




N 




N 




(K) 


(10'-' cm"^) 


(K) 


(10'= cm-^) 


EHVB 


58±4 


1.9±0.4 


160±18 


6.0±0.7 


EHV-R 


45±3 


2.4±0.5 


144±14 


7.6±0.8 


Broad 


46±2 


9.5±2.1 


82±4 


38.6±3.9 



3. Discussion 

3.1. Excitation conditions 

To determine the phy sical parameters, a smal l grid of RADEX 
slab models was run (Ivan der Tak et al ] l2007h . This grid is first 
used to determine CO excitation conditions by fitting r,ot (see 
Fig. IB. II in the on-line appendix). To reproduce Trot-^ISO K as 
observed in the EHV components, a kinetic temperature of 150 
K and density greater than -10* cm"^ is required. Alternatively, 
a kinetic temperature of ~500 K and density of ~10^ cm"^ can 
also reproduce the observed CO rotational temperature. In the 
following, both sets of values, (T, n) = (150 K, 10* cm"^; model 
1), and (500 K, 10^ cm"^; model 2), will be examined to con- 
strain the H2O excitation. Similar values are needed to repro- 
duce the lower CO rotational temperature of the broad compo- 
nent, r„t~80 K, (r, n) = (100 K, 10* cm-^) and (500 K, 5x10^* 
cm"^). The latter solution is excluded, however, because the en- 
velope density within the 20" beam is alw ays > 4x10^ £2Li- 
Inde pendent analysis of outflow emission dNisini et al. I [19991 
[2000) indicate that H2O is excited in a hotter, denser medium 
than CO, and for this reason a model with {T, n) - (500 K, 10* 
cm""*; model 3), is included in the analysis. 

To constrain the H2O excitation, the same physical condi- 
tions as determined for CO are taken as fixed and line widths of 
40 and 20 km s ' for the EHV-B and EHV-R components, re- 
spectively, and 50 km s ' for the broad component are adopted. 
C-type shock models show that CO a nd H2O cooling takes place 
at the same location in a shock ( Flowe r & Pineau des ForetsI 
[2OIO), thus validating this assumption. The H2O o/p ratio is as- 
sumed to be 3, and the only variable in the RADEX modelling 
is the total H2O column density. To constrain this, the ratio of 
the higher-excited H2*0 lines is taken with respect to the 2o2- 
lii line, and the calculated emission in this line is constrained 
to be greater than the observed value. The same is done for the 
Hj^O tentative detection and upper limit. A analysis is used 
to determine the best-fit value of the H2O column density. Once 



Fig. 3. H2O/CO ratio as a function of H2O column density for the 
three different models 1, 2 and 3. The best-fit results are marked 
with points. 

the column density is determined, the predicted intensity is com- 
pared to the observed intensity of the 2o2-l 1 1 line, and from this 
the filling factor is determined. The predicted column density is 
re-scaled to the average column density in the beam, and com- 
pared to N{CO) to determine the H2O/CO abundance ratio. 

The best-fit results are listed in Table |3] A^(H20) is found 
to be ~10'*-10'^ cm"^ for all components and all three mod- 
els, implying a beam filling factor of <0.01-0. 1 . This leads to an 
average H2O/CO abundance ratio of ~0.05-l .0 in all of the com- 
ponents; the range is caused by the excitation model rather than 
differences between the physical components. This result is very 
different from the rotational-diagram analysis above, where the 
abundance ratio is a factor of >100 lower The difference arises 
because H2O is not thermalised and optically thick, as opposed 
to CO, and thus non-LTE, radiative-transfer models are needed 
to treat the excitation properly. Figure |3] shows the H2O/CO ra- 
tio as function of A^(H20) demonstrating that only model 3 pro- 
duces abundance ratios <0. 1 . 

The RADEX models that match the H2O line emission 
best predict that the filling factor for the EHV component is 
very small, <0.01. Such low values agree with high spatial- 
resolution observations showing this component to be confined 
to the molecular jet and consisting of several sub-arcsec knots 
(e.g.. iHirano et al.H2010) . The model furthermore predicts that 
the ground-state transitions have optical depths in excess of 100 
(Fig. IB. 21 online appendix), and that the two Hj^O transitions 
and the higher-excited H2*0 transitions are also optically thick 
(t~3-4). Thus, the optical depth derived from the line ratios can 
only be considered a relative optical depth between the two tran- 
sitions. If the Hj^O lio-loi transition is assumed to be optically 
thin, a typical H2O column density of a few times lO'"" cm"^ 
is required. However, models with this column density predict 
that the intensity of the 2o2-lii line is close to that observed, or 
lower in the case of the broad component, implying filling fac- 
tors of unity inconsistent with the bullet nature. The CO emission 
remains optically thin (r<l) for all three models even with the 
small filling factors derived for H2O, consistent with the non- 
detection of '^^CO 3-2. 

3.2. Pliysical origin 

3.2.1 . EHV versus broad components 

Following the interpretation o f the EHV emission in IRAS04 191 
(ISantiago-Garcia et al.l 12009) it seems likely that the knots 
are internal working surfaces (IWS) along the jet that appear 
as strongly beam-diluted emission in single-dish observations. 
Regardless of the precise shock details, the high abundance of 
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Table 3. Modelled H2O line emission and inferred parameters. 





(lu cm ) 




Km s } 


— T— 

(xlO"-') 


HtO/CO 


Mnriel l'' ■ 

iVl*JLlCl I 












EHV-B 


2.4 


0.4 


3000 


1.3 


0.5 


EHB-R 


0.2 


0.4 


870 


4.4 


0.1 


Broad 


1.2 


1.5 


2070 


8.3 


0.3 


Model 2' : 












EHV-B 


4.0 


0.9 


2650 


1.5 


1.0 


EHB-R 


0.8 


0.4 


890 


4.3 


0.5 


Model 3'' : 












EHV-B 


0.8 


0.5 


3800 


1.1 


0.1 


EHB-R 


0.08 


0.9 


890 


4.3 


0.05 


Broad 


0.01 


0.2 


290 


65 


0.02 



Notes. Ortiio- + para-HjO. {T, n) = (150 K, 10'* cm^^) and (100 
K, 10* cm"^) for the EHV and broad components, respectively. (r, n) 
= (500 K, 10=^ cm-'). (r, n) = (500 K, 10** cm"') for all components. 



H2O in the bullets implies a high H2 fraction, since otherwise 
reactions with atomic H drive H2O back to OH and atomic O. 

The timescale for H2O to be produce d in gas as a function 
of temperature has been determined by B ergin et al.l ((1998,) for 
n=10^ cm"^. For 7=500 K, the timescale is ~100 years, which 
may be compared to the maximum dynamical time for a bul- 
let with velocity 50 km s"' to leave a beam with radius 11", 
240 years. Thus, it is possible to reform H2O from atomic gas 
over the timescales associated with the bullets. In the lower-T 
scenario, H2O will not form in the quantities observed here im- 
plying that the temperature is high in the bullets or that H2O is 
not formed in the gas phase. 

The gas-phase production of H2O is in agreement with 
model s of continuous spherical wind models of iGlassgold et alj 
(Il99lh . where the formation takes place from a dust-free atomic 
gas in sufficient quantities if the mass-loss rate exceeds 10"^ Mq 
yr"', which is to be compared to the time-averaged mass-loss 
rate in the jet of ~10"^ Mq yr"' (Hirano et al. 2010). The mass 
loss rate traced directly by the bullets is episodic (and thus higher 
than the time-averaged value) and the wind is not spherical, i.e., 
it is plausible that the H2O abundance is consistent with pro- 
duction in the IWS from atomic gas. To resolve the issue of the 
origin of the molecules and to break the degeneracy between the 
different models, it is necessary to obtain velocity-resolved data 
of chemically related species, such as O and OH. Furthermore, 
measuring the H/H2 ratio would test the models. 

The broad component, centered near ^lsr, is now ob- 
served in H2O emission in a number of low-mass sources 
dKristensen et al.l 12010. in prep.). It is interpreted as arising in 
the interaction between the outflow and envelope, and is seen so 
prominently in H2O emission (as opposed to, e.g., CO) due to 
the high production rate of H2O at this location. H2O is formed 
through both a series of neutral-neutral reactions pushing atomic 
oxygen into water, and through sputtering of icy grain mantles. 
Quantifying the contributions from the different mechanisms re- 
quire observations of other grain-surface products, e.g., CH3OH, 
and velocity -resolved observations of the intermediate gas-phase 
products O and OH. Tafall a et al.1 (I2OIOI) detect CH3OH emis- 
sion in both the EHV and broad components (but at an outflow 
position) indicating that a significant fraction of H2O is sputtered 
from grains. 



Table 4. H2O and CO cooling rates in 10"^ L, 







EHV-B 


EHV-R 


Broad 


Total 


H20: 


Obsei'ved emission 


0.4 


0.3 


1.1 


1.8 




Model 1 extrapol. 


11.0 


4.8 


26.6 


42.4 




Model 2 extrapol. 


10.7 


6.4 




17.1 




Model 3 extrapol. 


15.1 


6.5 


11.7 


33.3 


CO: 


Obsei'ved emission 


0.2 


0.2 


0.4 


0.8 




extrapol. 


1.1 


1.1 


1.8 


4.0 



3.2.2. Cooling 

The observed cooling rate in the three components is found by 
summing up the observed line emission (Table |4|l. For the EHV 
components, the H2O cooling rate is ~ 3.5x10 '* Lq, while it 
is ~1.3xlO Lq for the broad. Extrapolating the cooling to in- 
clude all lines with wavelengths greater than 60 yum gives a total 
H2O luminosity of 2^x10"^ Lq depending on which model is 
used (TablelU, i.e., the cooling by the lines observed with HIFI 
amounts to ~10%. The cooling is nearly the same for all three 
components demonstrating the power of spectrally resolving line 
profiles. This number may be compared to 4 .5x10'^ Lq as was 
found for H2O by ISO-LWS in an 80" beam dNisini et al. 200^), 
which applies to the sum of the components. 

The observed CO cooling rate is somewhat higher, ~2xl0"^ 
and 6.0x10"^ Lq, respectively. The CO cooling rate for the three 
components is extrapolated using the rotational diagrams to a to- 
tal of 4.0x10"^ Lq, lower than the H2O cooHng. The EHV com- 
ponents contribute ~50% to the CO cooli ng. ISO-LWS ob served 
a total CO cooHng rate of 3.0x10"^ Lq ('Nisini et al."2000), but 
in a much larger beam encompassing the extended outflow emis- 
sion and other physical components (Ivan Kempen et al.ll2010l) . 

4. Conclusions 

H2O emission readily traces the EHV gas that has previously 
only been seen in very deep integrations of other species, such as 
CO. The EHV components are attributed to shocks in the molec- 
ular jet, while the underlying broad component, now seen in sev- 
eral low-mass YSOs, is associated with the interaction of the out- 
flow and the envelope on larger scales. Thus, HIFI and H2O to- 
gether are ideal for revealing the dynamical components of low- 
mass star-forming regions. Both components are H2 rich, since 
otherwise the H2O abundance would be lower The two distinct 
components appear remarkably similar in terms of excitation and 
chemical conditions in spite of their different chemistries. They 
have equal contributions from H2O and CO cooling with H2O 
being the dominant coolant of the two. High spectral-resolution 
observations are a prerequisite for interpreting spectrally unre- 
solved PACS data of the same source. 
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Fig. A.l. (Top:) Gauss fit of the H2O 1 n-Ooo line (red) with the 
residual shown below (offset at -0.1 K). The central absorption 
feature is not fitted. (Bottom:) CO 3-2 scaled to the peak inten- 
sity of H2O in the bullets and overplotted on the H2O ln-Ooo 
Une. 

Appendix A: Gauss fit of data 

The data are fitted by Gaussians, with the exception of the nar- 
row absorption component. Figure lATI shows the example of the 
H2O 1 ii-Ooo line, with the residual plotted below. Furthermore, 
Fig. lA. II shows the CO 3-2 line scaled and overplotted on the 
H2O 1 1 i-Ooo line, highlighting the similar line profile in the EHV 
components. 

Appendix B: RADEX results 

The non-LTE, escape-probability code RADEX was used to cre- 
ate rotational diagrams for H2O and CO i n the optica l ly thin 
lim it. CoUisional rate coefficients are from iFaure et aP (l2007l) 
and lYang et al.l (l2010l) . respectively. The same transitions as ob- 
served were then used to calculate the rotational temperature as 
function of Tkin and n^. The results are shown in Fig. IB. II 

RADEX was also used to calculate the optical depth, t of 
the Hj^O lio-loi line at 557 GHz. Results are shown in Fig. 
IB. 21 where also the beam-filling factor is shown as a function of 
H2O column density for the three different models. 




10" 10= 10° 10' 



Fig. B.l. Tiot of optically thin H2O emission (top), optically thick 
H2O emission (A^^IO''' cm~^; middle) and optically thin CO 
emission (bottom) determined from RADEX simulations. T^ot is 
calculated from the same transitions as in Sect.|2] i.e., excluding 
the ground-state transitions. The gray area indicates the observed 
values of Trot. The points indicate the different physical condi- 
tions examined; model 1 (red) has {T, n) = (150 K, 10^ cm"^) 
for the EHV components and (100 K, 10^ cm"^) for the broad 
component. Model 2 (blue) has (T, n) = (500 K, 10^ cm^^) for 
the EHV components, respectively. Model 3 (green) has (T, n) 
- (500 K, 10^ cm"^) for all components. 




Fig. B.2. (Top:) t of the H^^O 1 lo-loi line. The black dashed line 
corresponds to the observed value in each of the three velocity 
components as specified in Sect. 2. (Bottom:) Beam-filling factor 
for the three different components and three different models as a 
function of total H2O column density. Best-fit results are marked 
with points in both plots. 
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Table 2. H2O, H2^0 and CO emission in each velocity component at the source position" . 



EHV-B EHV-R Broad 



Transition 


rms* 




j Tub ai) 


FWHM 






j Tub du 


FWHM 




^MB 


J Tmb do 


FWHM 


olsr 




(mK) 


(mK) 


(Kkms-') 


(kms-') 


(km s-') 


(mK) 


(Kkms-') 


(km s-') 


(kms-') 


(mK) 


(Kkms-') 


(kms-') 


(km s-') 


H2O 1 lo-loi 


9 


74 


2.54 


32.3 


-44.2 


122 


3.00 


23.0 


60.1 


294 


14.40 


46.0 


14.0 


1 1 ] -Ooo 


11 


101 


3.98 


37.1 


-41.1 


176 


4.00 


21.3 


60.1 


337 


14.76 


41.2 


20.1 


2o2-lll 


14 


93 


4.01 


40.3 


-43.7 


175 


3.86 


20.7 


59.8 


328 


17.59 


50.4 


18.0 


2ii-2o2 


14 


65 


2.47 


35.9 


-44.2 


109 


2.53 


21.9 


60.0 


207 


9.91 


44.9 


15.4 


3l2-3o3 


13 


50 


2.34 


44.0 


-45.4 


65 


1.14 


16.5 


60.7 


124 


7.76 


54.5 


16.3 


3 12-^21 


66 


147 


6.22 


39.7 


-43.1 


111 


3.32 


28.1 


55.4 


139 


7.94 


53.4 


14.8 




3 


4 


0.13 






5 


0.11 






5 


0.23 








11 




<0.25 








<0.19 








<0.28 






CO 2-1 


51 


184 


2.79 


14.2 


-51.6 


280 


4.12 


13.8 


60.3 


498 


38.53 


72.8 


0.8 


3-2 


39 


216 


4.82 


21.0 


-48.2 


467 


6.58 


13.2 


61.9 


657 


44.85 


64.1 


8.0 


4-3 


160 


467 


11.9 


24.0 


-45.1 


808 


15.1 


17.5 


65.7 


1023 


74.61 


68.5 


19.0 


10-9 


66 


269 


7.30 


25.5 


-47.2 


426 


8.42 


18.6 


61.3 


323 


17.77 


51.6 


10.1 



Notes. Obtained from Gaussian fits to each component. Measured in 1 km s ' bins. Obtained by fixing the position and width of Gaussian 
functions to the average value and only fitting the intensity. Upper limits are 3cr. 



Table 3, H2O, H^^O and CO transitions observed with Herschel-HIFI" . 



Transition 


V 


A 


EJkB 


A 


Beam 


'ml 


Hmb 




(GHz) 


(lim) 


(K) 


(10-3 s-') 


(") 


(min.) 




H2O lio-loi 


556.94 


538.29 


61.0 


3.46 


39 


13.0 


0.74 


lii-Ooo 


1113.34 


269.27 


53.4 


18.42 


19 


43.5 


0.75 


2o2-lii 


987.93 


303.46 


100.8 


5.84 


22 


23.3 


0.73 


2ii-2o2 


752.03 


398.64 


136.9 


7.06 


29 


18.4 


0.74 


3l2-3o3 


1097.37 


273.19 


249.4 


16.48 


20 


32.4 


0.75 


3 12-221 


1153.13 


259.98 


249.4 


2.63 


19 


13.0 


0.65 


H^'O lio-loi^ 


547.68 


547.39 


60.5 


3.59 


39 


64.3 


0.74 


lii-Ooo^ 


1101.70 


272.12 


52.9 


21.27 


20 


43.5 


0.75 


CO 10-9 


1151.99 


260.24 


304.1 


0.10 


19 


13.0 


0.65 



Notes. From the JPL database of molecular spectroscopy (iPickett et alJll998h . Total on + off + overheads integration time. Observed in 
the same setting as the main isotopologue. 



